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Abstract 



We consider a superprocess with coalescing Brownian spatial motion. We first prove a dual 
relationship between two systems of coalescing Brownian motions. In consequence we can express 
the Laplace functionals for the superprocess in terms of coalescing Brownian motions, which 
allows us to obtain some explicit results. We also point out several connections between such a 
superprocess and the Arratia flow. A more general model is discussed at the end of this paper. 
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In this paper we mainly consider the following branching-coalescing particle system which 
can be described intuitively as follows. A collection of particles with masses execute coalescing 
Brownian motions. In the mean while the masses for these particles evolve according to inde- 
pendent Feller's branching processes. Upon coalescing those particles involved merge into one 
particle with their respective masses added up. 

The above-mentioned particle system can be described using a measure-valued process Z. 
More precisely, the support of Z t represents the locations of those particles at time t, and the 
measure Zt assigns to each supporting point stands for the mass for the corresponding particle. 
This processes Z, referred as the superprocess with coalescing Brownian spatial motion (SCSM), 
was first introduced in [5]. It arises as a scaling limit of another measure- valued process, which 
was referred in [3] as the superprocess with dependent spatial motion (SDSM). As to SDSM, 
it arises as a high density limit of a critical branching particle system in which the motion of 
each particle is subjected to both an independent Brownian motion and a common white noise 
applied to all the particles. More precisely, the movement of the i th particle is governed by 
equation 



where (B^) is a collection of independent Brownian motions which is independent of the white 
noise W; see (3]. A similar model was also studied in |12j . 
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It was shown in Theorem 4.2 of [3] that, after appropriate time-space scaling, SDSM converges 
weakly to SCSM. A functional dual for SCSM was given in Theorem 3.4 of 5 . In addition, 
using coalescing Brownian motions and excursions for Feller's branching process, a construction 
of SDSM was found in [JJ], an idea that initially came from [I]. In this paper we always denote 
such a SCSM as Z . 

One of the most interesting problems in the study of a measure-valued process is to recover a 
certain dual relationship concerning the measure-valued process. Such a dual relationship often 
leads to the uniqueness of the measure- valued process; see ^U] for some classical examples on 
super Brownian motion and related processes. It is not hard to show the existence of Z as a 
high density limit of the branching-coalescing particle system. The main goal of this paper is to 
propose a new way of characterizing the measure- valued process Z via duality, in which the self 
duality for coalescing Brownian motions plays a key role. To this end, we first prove a rather 
general duality on two coalescing Brownian motions running in the opposite directions. We 
derive this duality from an analogous, essentially combinatorial, fact about coalescing simple 
random walk. With this duality we can express certain Laplace functionals for Z in terms of 
systems of coalescing Brownian motions. 

We could carry out some explicit computation thanks to the above-mentioned duality. In 
particular, we first show that, starting with a possibly diffuse initial finite measure Zq, Zt 
collapses into a discrete measure with a finite support as soon as t > 0. Then we can identify 
Zt interchangeably with a finite collection of spatially distributed particles with masses. When 
there is such a particle at a fixed location, we obtain the Laplace transform of its mass. The 
total number of particles in Zt decreases in t due to both branching and coalescing. When there 
is only one particle left at time t, we also recover the joint distribution of its location and its 
mass. Eventually, all the particles will die out. We further find the distribution of the location 
where the last particle disappears. Coincidentally, super Brownian motion shares the same near 
extinction behavior. 

Connections between superprocesses and stochastic flows have been noticed before. In |3] a 
superprocess was obtained from the empirical measure of a coalescing flow. Arratia flow serves 
as a fundamental example of coalescing flow. In this paper we point out several connections 
between Z and the Arratia flow. More precisely, the support of Zt at a fixed time t > can 
be identified with a Cox process whose intensity measure is determined by the Arratia flow. A 
version of Zt can be constructed using the Arratia flow. The general Laplace functional for Z 
can also be expressed in terms of the Arratia flow. 

Replacing the Feller's branching process by the square of Bessel process, we discuss a more 
general model at the end of this paper. The mass-dimension evolution of such a model can also 
be characterized by coalescing Brownian motions. 

The rest of this paper is arranged as follows. As a preliminary, we first state and prove a dual 
relationship on coalescing Brownian motions in Section [2j In Section 12 we define the process Z 
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as a weak limit of the empirical measure for the branching-coalescing particle system. Then we 
proceed to prove the duality between Z and coalescing Brownian motions. The uniqueness of Z 
follows from such a duality immediately. We continue to study several properties of this process 
in Section^ We further discuss the connections between the Arratia flow and Z in Section 03 
At the end, we propose a more general model and establish its duality in Section 6. 

2. Coalescing Brownian motions and their duality 

An m-dimensional coalescing Brownian motion can be described as follows. Consider a system 
of m indexed particles with locations in R that evolve as follows. Each particle moves according 
to an independent standard Brownian motion on R until two particles are at the same location. 
At this moment a coalescence occurs and the particle of higher index starts to move together 
with the particle of lower index. We say the particle with higher index is attached to the particle 
with lower index, which is still free. The particle system then continues its evolution in the same 
fashion. Note that indices are not essential here, the collection of locations of the particles is 
Markovian in its own right, but it will be convenient to think of the process as taking values 
in R m rather than subsets of R with at most m elements. For definiteness, throughout this 
section we will further assume that the particles are indexed in increasing order of their initial 
positions: it it clear that the dynamics preserve this ordering. Call the resulting Markov process 
X = (X\, . . . , X m ). 

Write 1{B}(.) for the indicator function of a set B. The distribution of X(t) is uniquely 
specified by knowing for each choice of y\ < y 2 < ■ ■ ■ < y n the joint probabilities of which 
"balls" X 1 (t),X 2 (t), . . .,X m (t) lie in which of the "boxes" ]y x , 2/2], ]U2, Va\, ■ ■ ■ , }y n -i, Vn]- That 
is, the distribution of X(i) is determined by the joint distribution of the indicators 

i£(t,y) -=HMt) e]i/j,v;+i]} 

for 1 < i < m, 1 < j < n - 1 and y = (y 1 , . . . , y n ). 

Suppose now that Y := (Y\, . . . , Y n ) is another coalescing Brownian motion. The distribution 
of Y(i) is uniquely specified by knowing for each choice of x\ < x 2 < ■ ■ ■ < x n the distribution 
of the indicators 

J£"(t,x) := l{ Xl e]Yj(t),3S +1 (t)]} 

for 1 < i < m, 1 < j < n — 1 and x = (x±, . . . ,x m ). 

The next "balls-in-boxes" duality is crucial in characterizing the distribution of the measure- 
valued process concerned in this paper. 

Theorem 2.1. Suppose in the notation above that X = (X±, . . . ,X m ) is an m- dimensional 
coalescing Brownian motion and Y = (Y±, . . . ,Y n ) is an n-dimensional coalescing Brownian 
motion. Then for each t > the joint distribution of the m x (n — 1)- dimensional random array 
(I^(t, Y(0))) coincides with that of the m x (n — 1)- dimensional random array X(0))). 
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Theorem l2.1| generalizes Theorem 1.1 in |16j . Some other more elaborate dualities on coalescing- 
reflecting Brownian systems can be found in |14j and |13j . 

We first prove the counterpart of Theorem 12.11 for continuous time simple coalescing random 
walks, which is interesting in its own right. Notice that X is a coalescing Brownian motion 
if and only if X,- t is a (.F x )-Brownian motion for each 1 < i < m, and (Xj — Xi)/y2 is a 
(J^ x )-Brownian motion stopped at 0, where (J^ x ) denotes the filtration generated by X. Then 
Theorem 12.11 follows from a straight forward martingale argument proof of the convergence of 
scaled random walk to Brownian motion. 

A p-simple random walk on Z is a continuous time simple random walk that makes jumps 
at unit rate, and when it makes a jump from some site it jumps to the right neighbor with 
probability p and to the left neighbor with probability 1 — p. An m-dimensional p-simple 
coalescing random walk is defined in the same way as the coalescing Brownian motion at the 
beginning of this section. When p = ^ we just call this particle system a simple coalescing 
random walk. 

Some notation is useful to keep track of the interactions among the particles in the coalescing 
system. Let V m denote the set of interval partitions of the totality of indices N m := {1, . . . , to}. 
That is, an element tt of V m is a collection tt = {Ai(ir), . . . ,Ah(ir)} of disjoint subsets of N m 
such that Ui^C 71 ") = and a < b for all o G Ai, b G Aj, i < j. The sets Ai(ir), . . . Ah(ir) 
consisting of consecutive indices are the intervals of the partition tt. The integer h is the length 
of 7r and is denoted by 1(tt). Equivalently, we can think of Vm as a set of equivalence relations 
on N m and write i ~„- j if i and j belong to the same interval of n G V m - Of course, if i ^ j, 
then i ^ k ^ j for all i < k < j. 

Given tt E V m , define 

ai(ir) := min^4j(-7r) 
to be the left-hand end-point of the i th interval Ai(ir). Put 

:= {(xi, . . . , x m ) G Z m : x\ < . . . < x m and X{ = Xj if i ~„- j} 

and 

Z™ := {(xi, . . . , Xm) £ Z m : x\< . . . < x m and Xi = Xj if and only if i ~ 7r j}. 

Note that Z m is the disjoint union of the sets Z^, tt G V m . 

Write X = (X%, . . . , X m ) for the coalescing random walk. If X(t) G Z™, then the free particles 
at time t have indices ai(7r), . . . , au^\{Tr) and the i th particle at time t is attached to the free 
particle with index 

min{j : 1 < j < m, j ~ w i} = max{a fc (7r) : a k (ir) < i}. 

In order to write down the generator of X, we require a final piece of notation. Let {e^ : 1 < 
i < k} be the set of coordinate vectors in Z fc ; that is, is the vector that has the i th coordinate 
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1 and all the other coordinates 0. For tt G V m , define a map 

K n (x) = K 7T (xi, . . . , X m ) := ^ ai (7r), • • • , x ai^(ir)j 

Notice that K n is a bijection between Z™ and {x G Z'W : Xl < £2 < • • • < x^}, and we write 
K~ l for the inverse of K n . For brevity, we will sometimes write x^ for K n (x.). 

Write B(Z m ) for the collection of all bounded functions on Z m . The generator G of X is the 
operator G : B(Z m ) -> 5(Z m ) given by 

G/(x) := p £ / o K-^x, + ef >) + (1 - p) £ / o K~\^ - e' W ) 
i=i i=i 

- /(tt)/ o tf" 1 ^), / G £?(Z m ), xd™, f£ 7> m . 

This expression is well-defined, because if x € Z™, then x^, x^ + e'^ and x T — e'^ are all in 
{x G Z'W : xi < x 2 < ... < 

Note: From now on we will suppress the dependence on dimension and write e-^ as e^. 

Write Z' := Z + ^ = {i + 5 : j £ Z}. An n-dimensional g-simple coalescing random walk on 
Z /n and its generator H can be defined in the obvious way. Such a process, with q = 1 — p, will 
serve as the process dual to the p-simple coalescing random walk on Z m in the following way. 

Fix x G Z m with x± < . . . < x m and y G Z /n with y 1 < . . . < y n . Put 

I£(t,y) :=l{Xi(t) €] Vj ,y j+1 ]} 

and 

/£"(<, x) := e]Y 3 (t),Y J+1 (t)]} 
for 1 < i < m and 1 < j < n — 1. 

Lemma 2.2. Suppose in the notation above that X = (X±, . . . ,X m ) is an m- dimensional Z m - 
valued p-simple coalescing random walk and Y = (Yi, . . . ,Y n ) is an n-dimensional U n -valued 
(l—p)-simple coalescing random walk. Then for each t > the joint distribution of the mx (n — 
1)- dimensional random array (i^*(i, Y(0))) coincides with that of the m x (n — 1)- dimensional 
random array (J^"(t, X(0))). 

Proof. For a function g : {0, l}" 1 *™- 1 ) _ > R ) a vector x G Z m with x\ < . . . < x m , and a vector 
y G Z /n with yi < . . . < y n , set 

y) := 9 (i{]yi, y2]}(zi), • • • , i{]y n -i, yn]}(zi), • • • , i{]yi, 2/2]} (zm), . . . , i{]y„_i, y„]}(x m )) . 

We may assume that X and Y are defined on the same probability space (fi, J 7 , P). We need to 
show that 



(2.1) 



P[S(X(t);Y(0))]=P[S(X(0);Y(t))]. 
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For x G Z m , put <7x( - ) := #( x ;')) an d f° r Y e ^'™> P u ^ 5y(') := ff("?y)- in order to establish 
(|2.1I) . it suffices by a standard argument (cf. Section 4.4 in [:6J to show that 

(2-2) G(5y)(x)=F(5x)(y) 

(recall that G and i7 are the generators of X and Y, respectively). 
Fix x G Z™ and y G Z'™ for some vr G P m and w G TV Put 

/+ : = {f : 1 < i < /(vr), x ai(7r) + - = y a ^ m) for some 1 < j < l{m)} 

and 

1 

I :={i:l<i< /(vr), x ai(n) - - = y a ^ m) for some 1 < j < l(w)\. 

Similarly, put 

J~ ■= {j ■ 1 < j < 'H, - - = x^ (7r) for some 1 < i < /(vr)} 

and 

J+ : = {j : 1 < j < /( CT ), y aj(ro) + - = x Qi(7r) for some 1 < i < /(vr)}. 

Recall that z Ql(7r) < . . . < a^,^) and V ai {vj) < • • • < 2/a i(ro) ( CT )- Therefore, for each i G 7+ 
there is a unique j G J~ such that x^.^) + \ = y a] ( m ) and mce versa. Fix such a pair (i, j). 
Observe that 

x':=x+ *k=K-\^ + e i ) 

and 

y':=y- £ e£ = ^(y ro - e,-). 

We are going to verify that 

(2.3) (i{]^,^+i]}(4)) = {HWfiy'j'+iftM) 

by considering all the possible scenarios. 
Given any i' G Ai(Tr) we have: 

• for f = ctj(w) — 1, 

l{]yy, Vf+i]}M = !{]%', Vj>+l]}(xi> + 1) 
= 

= Vi'+i - i]}(zi') = i{]yi' 5 yi'+i]}(^')> 

• for aj{m) < f < maxAj(w), 

yj'+i]}(4) = Hlyi'^j'+i]}^' + !) 

= o 

= !{]%■' - l.yj'+i - i]}(af»0 = H}y'j',yj'+i}}(xi'), 



• for j' = max Aj(w), 

HW^Vj'+i]}^') = yf+i\}( x i' + !) 

= 1 

= l {\vj' - hyj'+i]}(xi>) = HWj>,y'j'+i}}(xi>), 

• and for j' < ctjiva) — 1 or j' > max Aj {w) , 

= 

= !{]%•', yj'+i]}(^i') = 

Moreover, given any z' Ai(n), we have a;^ / x a .r^y Hence 

• for j' = atj(m) — 1, 

H]yf,yj'+i]}(xi>) = l{\Vj'>yj'+i]}M = lQyj',y?+i - l]}(av) = 

• for j' = m&xAj(w), 

• for aj(ro) < / < max Aj(w), 

!{]%•'» %'+i]} (40 = %-'+i]}(«t') = !{]%•' - - = l{]l^yJ-' +1 ]}(2v), 

• and for j' < a.j(zu) — 1 or j' > max Aj(ro), 

Combining the above observations yields (|2.3j) . 
Therefore, 

9 y ° ^Q^O^r + e») = g x o K" 1 ^ - a,-). 
Furthermore, it is easy to see for i' / + that 

9y ° ^^Xtt + ej/) = <? y o iir~ 1 (x 7r ) 

and for / that 

Similarly, for any i G /" there exists a unique j G J + such that — \ = y aj (m) arm 

versa. For such a pair we have 

<?y ° ^tT 1 ^ - e») = 5x o K^^yro + e,). 
Furthermore, we see for i' I~ that 

5y ° ^tT 1 ^ - ej/) = 5y ° ^Hxtt) 

and for j' J + that 

<?x ° if^Cy^ + e,') = g x o ET~ 1 (y ro ). 
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Lastly, note that 

9y ° K~ x {-x^) = <?(x;y) = g x o K m l {y m ) 

and so 

G(g y )(x)-H(g x )(y) 

= P^2{9y° K~ x (x n + ei) - g y o jf" 1 ^)) 

i=l 

+ (! " P) S ° - ei ) - o IT" 1 for)) 

i=l 

- P (# x ° ^Hy^ - ej) -g x o K~ 1 {y w )) 
3=1 

-(l-p)J2 (g x o K~\y w + ei) - g x o i^y-)) 
J=l 

= P^5y° -^(xtt + ei) - p ^2 5x ° A'-^y^ - ej) 
iei+ jeJ- 

+ (1 " P) Yl 9y ° ^'(x. - ^) - (1 - p) 5x o K-\y w + e,) 
= 0, 

as required. □ 

Remark 2.3. For discrete time coalescing simple random walks such a duality is evident from 
Fig. 7 in 13 j. But the duality seems to be less apparent for continuous time coalescing simple 
random walk. 

3. Existence and uniqueness 

A construction of Z was given in jS] using Feller's branching excursions. In this paper we 
adopt a weak convergence approach, which is commonly used in the study of measure-valued 
processes. 

Recall that a nonnegative valued process £ is a Feller's branching process with initial value 
x > if it is the unique strong solution to the following stochastic differential equation 

,-t 



Jo 



where 7 is a positive constant and B is a one-dimensional Brownian motion. is a martingale. 
It has a Laplace transform 

2\x 



(3.1) P[exp{-A£(t)}] = exp 



2 + X-ft 
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its extinction probability is given by 

2x 



P{£(t) = 0} = exp 



jt 



see Section II. 1 and II. 5 in |1U| . 

Observe that independent Feller's branching processes are additive; i.e. if £ and rj are two 
independent Feller's branching processes (with the same parameter 7), then £+77 is also a Feller's 
branching process. This fact will be used repeatedly in our discussions. 

Write Mp(R) for the space of finite measures on R equipped with the topology of weak 
convergence. Given any finite measure Z$ on R, put z := Zq(K). For any positive inte- 
ger m, let (£1 , •••,£m ) be a collection of m independent Feller's branching processes each 
with initial value z/m. Choose {x\ , . . . , x n ) to be i.i.d. samples from distribution Zq/z. Let 
(x[ m \ . . . , X^ ) be an m-dimensional coalescing Brownian motion starting at (x\, . . . ,x m ). 
Moreover, we always assume that (£,[ m \ . . . ,£^T' ) ) and (x[ m \ . . . ,xffl) are independent. Let 
5 X denote the point mass at x G R. Then 



x\ m >{t) 
i=i 

defines a M^(R)-valued process. From now on we will suppress the dependence of m in £ ? - 
and^ (m) . 

Recall that a collection of processes {Z a ,a G 1} with sample paths in Z)(Mf(R)) is C- 
relatively compact if it is relatively compact and all its weak limits are a.s. continuous. The 
proof of the next lemma is standard; see, e.g. the proofs for Lemma 3.2 in ^S] and Proposition 
II.4.2 in HOI. 

Lemma 3.1. {Z^} is C-relatively compact. 

Proof. We first check the compact containment condition. For any e > and T > 0, choose a 
compact set K C D(R) such that F{X X G K§} < e 2 . Let K := {x t : x G K ,t < T}. Then K 
is compact in R, and 

P{Xi(i) G K c , 3t < T} < P{Xi G K%\ < e 2 . 

Write 

iV m := #Ik :=#{!< i< m: X^t) G K°, 3t < T}, 

where #Ik denotes the cardinality of the indices set Ik- Conditioning on N m , by the additivity 
for Feller's branching processes, we see that Y2i<=i K is a Feller's branching Process with initial 
value N m z/m. Then by Doob's maximal inequality, 



sup > y £i(tj > e 

0<t<T 



N m z 



me 



10 

Therefore, 



sup z[ m) {K c ) > e \ < F{ 

0<t<T 



sup 

0<t<T 



E > e U 



P[iV m ]z 



me 



For any / € C^(R) put zj: (f) := f{x)z[ {dx). Now we are going to show that 
{Z( m )(_/)} is C-relatively compact in D(M). By Ito's formula, we have 



4 m) (/) = E 



i=i 



m 



+- 



The additivity for gives 



sup 

0<s<t 



i=l 



< oo and P 



sup E &( s )£i( s ) 



0<s<t 



*J=1 



< oo, t > 0, 



then X/^=i Jo £i( s )f"(Xi(s))ds is C-relatively compact following from the Arzela-Ascoli theorem 
and Proposition VI. 3. 26 of [7J. 
Note that 

/ m \ m -t 

(E/ eiW/ww^w = E / ei(^(s)/ , (x i ( s ))/ / (x i ( s ))d(x i ,^) SJ 

where {Xi,Xj) s = s — A s and Ty := infjs > : Xi(s) = Xj(s)}. By Arzela-Ascoli theorem 
again, {{J2^L\ fo£,i( s )f'(Xi(s))dXi(s)).} is C-relatively compact. Theorem VI. 4. 13 and Proposi- 
tion VI. 3. 26 in [JJ then imply that the collection of martingales {Y^iLi Jo £i( s )f'{Xi(s))dXi(s)} 
is C-relatively compact. 

Similarly, {YliLi In f(Xi(s))d£i(s)} is also C-relatively compact. Moreover, 



- m 

— /~2f( x 'i) ~> z o(f) a.s.. 



{Z^ m ^(/)} is thus C-relatively compact. Consequently, by Theorem II. 4.1 in ^U] we can conclude 
that {Z^} is C-relatively compact. 

□ 

Write Z for the weak limit of {Z^}. The Laplace functional of Z can be obtained from the 
duality in Theorem 12.11 As a result, its uniqueness is settled. 

In the sequel we always write (Y\,...,Y2 n ) for a coalescing Brownian motion starting at 
(yi, V2n) with yi< ... < y 2n - 
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Theorem 3.2. {Z^} has a unique weak limit Z in C(Mp(R)). Given aj > 0, j = 1, . . . ,n, 
for any y\ < 1/2 < • • • < yi n and any t > 0, we have 



(3.2) 



exp <{ -^a J Z t (]y 2 j-i,y2j]) 



cxp 



Z (dx) 



2E?=iOjl{]^-iW.^i(t)]}(a 



Proof. We first condition on (£i(t), . . . , £ m (t)). By Theorem 12. II we can show that 



(3.3) 



j=l i=l 



ex p i -EE^w^^hW-^W]}^ 

i=l i=l 



(6(0) 



where (Yi, . . . , y^n) is independent of (Xi) and (£,). 

Now take expectations on both sides of (|3.3j) and then condition on (ccj) and (y(i)). Since 
. . . ,£m(t) are independent of each other, and they are independent of (xi) and (Yi(t)), it 
follows from Q3.1j) that 



exp { ~ Yl a j Z t m) (]y2j-i,y2j}) 



(3.4) 



ex P { -zZT, a Mt)n]Y23-i(t),Y 2j (t)]}(x t ) 
3=1 i=l 



n 



cxp 



f 2zE" = i"ii{]^-i(t),y2i(t)]}(x,) 
[ m{2 + 7iE?=i%i{]y^-i(<),y^(<)]}(x0} 
, f 2zE^i^i{]y 2 ^iW,y 2j (t)]}(x) 

Zq{cIx) exp 



m(2 + 7*Ei=i oii{]y%--i(t) s y^(*)]}(x)) 

Let Z be any weak limit of {Z^}. Let m — > 00 in (|3.4jl . Then 



lim P 

m— >oo 



exp { ~ Yl a j z t m) (]y2j-i,y2j]) 
3=1 



(3.5) 



lim P 

m— >oo 



cxp 



Zo(dx 

) 

Zo(dx 



2E; =1 « J i{]y 2j -i(t),y 2j (t)]}(x) 

j m(2 + 7* E"=i a J -l{]y ?J -_ 1 (t), y 2j (i)]}(*)) 

2E-=i«ii{]y2i-i(i),y2iW]}(^) 



2+7*ELi%i{]y2i-i(*),y2i(t)]}(^) 
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We still need to make sure that 

IF" exp l-Y^aj z t(]V2j-i,y2j]) 



i 


= lim P 


exp | 




m— >oo 





Am), 



To this end we can suppose that yi < y 2 < • • • < 2/2n- Then for small enough e > 0, similar to 
1)3. 4 J) we have 



(3.7) 
exp ^ 

< 1 - 







✓ 


I] 


-p 


exp < 

v. 



- X] a i z i m) (]y2j-i + e, y 2 j - e]) 

3=1 

expj-j: ajZ t (m) (]y 2 j-i - e, 2/2j-l + e]U]y 2 j - e, y 2 j + e]) 



^ajZ t (m) (]?/ 2i _i - e,y 2i + e]) 
i=i 



1 

-,2n t^ri 



Zo(dx) exp 



-2z-E- =1 a J l{]^- 1 (t),^_ 1 W]U]^.(t) ; y 2 ^(t)]}(x) 
m(2 + 7< E;=i « 3 1{]^-i W, W]U]^ W, ^(t)]}(x)) 



<i-p{n^ 1 {y/(i) = y/(i)}} 1 

where {Y{, Y", ■ ■ ■ , ^2n' ^2n) ^ s a coalescing Brownian motion starting at (yi — e, y\ + e, . . . , y 2n — 
e>2/2n + e)- Clearly, 1)3 .7|) converges (uniformly in m) to as e — » 0+. So, 1)3 .6|) holds. 

It is clear that the distribution of Z is uniquely determined by 1)3.2)1 . So, Z is the unique weak 
limit of {Z^}. 

□ 

The moments of Z can be obtained immediately from ()3.2)) . 
Proposition 3.3. For any yi < y 2 < . . . < y 2n and t > 0, we aawe 



and 



z t ( X^fej-i^]} 



z 2 \ Y,^{]Y2i-i{t),Y 2j m 

i=i 



Z (^a J l{]F 2 ,_ 1 (t),y 2i (t)]} 



+ 



i=i 



Remark 3.4. Another consequence of duality 1)3.2)1 is that Z is a Markov process; see Theorem 
3.3 in ^B] for a proof on a similar model. 
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Martingale problem is often used to characterize a superprocess. Z is the solution to the 
martingale problem (see [5]): for any <fi S C 2 (IR), 

1 /"* 

M t {<j>) = Z t {4>) - Z (<t>) - - / Z s (<jJ')ds, t > 0, 

z Jo 

is a continuous martingale relative to {Ft)t>o with quadratic variation process 

{M{4>)) t = 1 f ' Z s (<p 2 )ds+ f als [ 4>\x)^{y)Z s (dx)Z s {dy), 
Jo Jo Ja 

where A = {(x, x) : x G M}. 

But a remarkable feature of such a martingale problem is that its solution is not unique. 
For example, let £i and £2 be two independent branching processes each with initial value 1. 
Let B\ and B2 be two independent Brownian motions. Assume that (^1,^2) and (-61,-82) are 
independent. Then Z' t := Ci(*)^Si(t) +^2{t)^Bi(t) i s another solution to this martingale problem; 
also see |16j for a similar counter example. 



4. Some properties 
Our first result in this section is a straight forward consequence of Theorem 13.21 

Proposition 4.1. For any y\ < yi < . . . < y2 n and t > 0, we have 



(4.1) 



exp < 



-\J2 z t(]y2j-i,y2j]) 



2A 



cxp 



2 + A 7 t ^ 

.7=1 



Y,Zo(]Y2j-i(t),Y 2l (t)]) 



A > 0. 



Proof. Observe that YlJ=i l{]^2j-i(0> Y2j(t)]}(x) is either or 1, then (|4.1|) follows readily from 
(IT21 . 



□ 



Proposition 14.11 allows us to carry out some explicit computations on Z. First, by letting 
A — > 00 we can easily see that 

2z 



P[Z t ((-oo,oo)) = 0] =exp|- — 

We are going to further study the probability that Zt does not charge on an arbitrary finite 
interval. For any x,y,a and b, put 



x - y _ 
x := — —,y 



V2 



x + y a + b ~ b — a 

, a := — =- and b :- 



V2 



V2 



V2 ' 
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Proposition 4.2. Given a < b and t > 0, we have 

1 

2vrt 



P{Z t (]o,6]) = 0} 



dx 



o 



(4.2) 



exp 



dy exp 

(y-b) 2 



2Z (]x,y]) {x-af 



tj 2t 
< (y + b) 2 \ 



2t 



exp 



2f 



+ 



'2-Kt Ji 



dx exp 



x 
2t 



Proof. Let A — > oo in (j4,l[) we have 



'{Z t (]a,b]) = 0} 



exp 



exp 



—zoimt)^)}) 



it 



Z (]Y 1 (t),Y 2 (t)])y i Y 1 (t)^Y 2 (t) 



+ F{Y 1 (t) = Y 2 (t)}, 



where (Yi,Y 2 ) is a coalescing Brownian motion starting from (a,b). 

To find the distribution of (Yi, Y2), one could rotate the coordinate system anti-clockwise by 
7r/4. Under the new coordinate system (Yi,Y 2 ) becomes a process (Yl^Y^) such that Y[ is a 
Brownian motion starting at 5, 1^' is a Brownian motion starting at b and stopped at 0, and Y{ 
and are independent. So (|4.2|) just follows from the reflection principle for Brownian motion. 

□ 

Write St for the support of Zf. Intuitively, starting with particles of a total initial mass 
Zq(R), as soon as t > the particles near —00 and 00 will die out due to branching. Z t is then 
expected to be supported by a finite set because of coalescence. The next two results concern 
the cardinality of St- 

Proposition 4.3. Given a < b and t > 0, we have 
P[#S t n]a,b}} 



(4.3) 



1 



dz 



dx / dyyexp 



2Z Qx, 



V2zf 



1 id V2nt 2 

Proof. It is easy to see from (|4.2j) that for any 2£l, 
F{Zt(dz) ± 0} 



21: 



(4.4) 



dz 



dz 



7Tt V2~1Tt 2 



dx 



dyy exp 



2Z (]x,y}) (x-V2z) 2 + y 2 



fry 



2t 



Then (|4.3|) is obtained by taking integrations on both sides of (|4.4jl from a to b. 



□ 



Proposition 4.4. With probability 1, #St < 00, Vi > 0. 
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Proof. Given s > 0, we first claim that F[#S S ] < oo if Z has a bounded support. Suppose that 
Z (] - b, b}) = 1 for some b > 0. Then by g3J), 



y2ixs J-oo J-oo Jo \ I s 7 J / [ 2 « 

1 r^r^.^f H*^n r_£ 



IlISS J-oo Jo V I s 7 J / I 2s 

J /-oo ry+V2b f y 2 



< —==- I dy dxy exp < 

\/2lTSS Jo J-y-V2b [ 2s 



< OO. 

Our claim is proved. 

Now given any integer j, let Tfj(s) be the Feller's branching process with initial value 77,- (0) := 
Zo(]j,j + !])• Since 

EPWW-t( 1 -o,{-»}) S f»-?, 

jr' = — OO j = — O0 j = — OQ 

by Borel-Cantelli lemma we have that, with probability 1, ryj(s) 7^ for only finitely many values 
of j. 

Therefore, for any i > 0, with probability 1, Z t / 2 must have a bounded support. The Markov 
property for Z, together with the claim from the first part of the proof, implies that #St < 00 
a.s.. 

Finally, by the Markov property for Z we conclude that P{#St <oo,Vi>0} = l. □ 

By Proposition 14.41 as soon as t > 0, St becomes a finite set. For any z G St, we associate it 
with a particle located at z with mass Z t ({z}). We can thus identify Z t interchangeably with 
a collection of spatially distributed particles with masses. As time goes on, the total number 
of particles decreases either because two "alive" particles coalesce into one particle, or because 
each particle disappears due to its branching. 

Since #St < 00, a small neighborhood of z contains at most one particle in Zt- When there 
is such a particle, we want to find the distribution of its mass. Formally, we are looking for 

F[exp{-XZ t ({z})};Zt({z}) > 0] . 

Proposition 4.5. For any zeR and t > 0, we have 
F[exp{-XZ t (dz)};Z t (dz) > 0] 



(4.5) v^Tt 2 7-oc Jo V I 2 + A 7 t J r \ fry 

exp 



(x - V2z) 2 + y 2 



2t 
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Proof. We fix (&(*)) first. Apply Theorem |2~T1 to 

( m ~\ ( m 

exp -A^ei(*)l{]o,6]}(^(*)) 1 5>(t)l{]M]}(*<(t)) > 



. i=l 



Then condition on (Yi(t), Y^i)) and take an expectation with respect to (£»(*))• Similar to the 
proof for Theorem lcS.2l we have that 



(4.6) 



i=l J i=l 

( m 

exp\-\Y,mH)Yi(t),Y 2 (t)]}(; 



8=1 



1=1 



^^(i)i{]y 1 (i),y 2 (i)]}(x i ) = o 



,i=i 



2A4 m) (]y 1 (^),r 2 (Q]) 

2 + Xjt 



cxp 



2Z^(]F 1 (t) ; F 2 (t)]) 
7 t 



exp 



Therefore, 



F[exp{-XZ t Qa,b])};Z t Qa,b]) > 0] 

2AZ (]y 1 (t),r 2 (t)]) | 



exp 



(4.7) 



2 + A 7 t 



cxp 



(x — a) 2 
2* 



exp 



exp 

2 AZ (]£,£]) 
2 + A 7 t 

(y - ft) 2 1 

2t | 



2Z (]y 1 (t),y 2 (t)]) | 
it J 

2^o(]5,jf]) 



cxp 



cxp 



fry 

(y + fr) 2 ) 

2t I 



So, (|4.5f> is obtained by letting 6 — ► o+. 



□ 



At a fixed time t > 0, with a positive probability there can be only one particle (with a 
positive mass) left. When this happens, we are interested in the joint distribution of the mass 
and the location of that particle. More precisely, we want to find 

F [exp{-AZ t (M)} ; Z t (R) ^0,S t C dz] . 

Proposition 4.6. For any z£l and t > 0, we have 
F [exp{-AZi(R)} ; Z t (R) ^0,S t C dz] 



dz 



(4.8) 



V27rf 2 

dz 



dx / dyyexp 



2XZ (]x,y]) 2Z (]x,y] c ) (x-V2z) 2 + y 2 



2 + X^t 



jt 



2t 



2z 



Proof. Put 



B := I J2ti(t)H}Yi(t),Y 2 (t)} c }(xi) = \ for Xi := X^O). 



It follows from Theorem 12, II that 
(4.9) 



exp ! -\J2Ut) [ ;X>(t)l{]a, &]}(*»(*)) > ? X;^(t)l{]a,&] c }(X,(t)) = 

I i=l J i=l i=l 

expl-\J2Ut) \ -,^i(t)H]Yi(t),Y 2 (t)]}(xi) >0,B 

I i=l J i=l 

( m ^ m 

exp \ -xj2mn}Yi(t),Y 2 (t)}}(x l ) k-^CiWi{]n(t),y 2 (t)]}(^) > 0,5 



i=i 

m 



i=l 



exp ^ - A £ Ci (*) ,Y 2 (t)} }( Xi ) \; B 



i=i 



£e i (t)l{]yi(t),r 2 (t)]}( 3 ; i ) = 0;B 



i=l 



exp 



2AfEI=il{]^iW^2(i)]}(^)l „„„ f 22E£il{]n(*),l r 2 (*)] c }(^) 



m(2 + A7t) 



exp 



rwyt 



i=l 



where in obtaining the last equation we have used the fact that, given (Yi(£), 5^(£)) an d (a^), 
expJ-A^Ci(*)l{m(t)T2(t)]}(^) [ and event i ^(^{^(t), Ya^)]^;) = I 



i=l 



. i=i 



are independent. 

Now letting m — > oo in (|4.9|) we have 

F [exp{-AZ t (M)} ; Z t (R) ^0,S t C (a, &)] 

2XZ (]Y 1 (t), Y 2 (t)}) 2Z (]Y 1 (t), Y 2 (t)) c ) 
2 + Xyt -ft 

2XZ (]x,y}) 2Z (]x,y]' 



exp 



exp 



(4.10) 



1 

2vrt 



oo poo 

dx \ dyexp 

-oo JO 



exp 



H — p= exp 



dx exp 

Finally, (|4~H|) is obtained by letting b —> a+ in (|4.1U|) 



2t 

it] h 



2 + A7i 7t 



exp 



2z 
~~t 
(x — a) 2 
2t 



exp 



2t 



exp 



2z 



Remark 4.7. Let A = in (|4.1U|) . We then obtain a result on the range of St- 
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The total number of particles in Z will decrease one by one. Put 

r := inf{s > : #S S = 1}. 

Then r < oo is the first time when there is exactly one particle left. The distribution of r is 
given in the following Proposition. 

Proposition 4.8. 



. {T<t}=/ 4 i r^r iw ^ 2z ^> I'-w+v 1 



— oo 



1 f 2zl\ , f 2z 

exp < > + 1 - exp 



(4.11) 

Proof. Observe that 

/oo 
P{Z f (M) ^0,S t Cdz} + F{Z t (R) = 0}, 
-oo 

then (|4. 1 1|) follows from Proposition 14.61 
Let 

T : = inf{t > : Z t (R) = 0}. 
T is the time when all the particles disappear. Its distribution can be found easily. 

P{T <t} = ¥{Z t (R) = 0} = exp |-^| . 

Let F denote the location of the last particle immediately before extinction, i.e. {F} = St- 
We could recover the explicit distribution for F. 



□ 



Proposition 4.9. F has the same distribution as Xt, where X is a Brownian motion with 
initial distribution Zq/z, and X and T are independent. 

Proof. First assume that 

m 

Zo = ^2 ai5 Xi 

i=l 

with dj > and Y^!i=i a i = a - Then 

m 

z t ■= ^2£i(t)8xi(t)> 
i=i 

where £i(0) = dj > and (X\, . . . , X m ) is a coalescing Brownian motion starting at (x\, . . ■ , x m ). 
Write Ti := inf{i > : £i{t) = 0}, i = 1, . . . , m. Then T := maxi<j< m Tj. Therefore, 

m m 
i=l i=l 
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Our first observation is that 

P{T, <t}= P{&(i) = 0} = exp 

Then P{T = T{\ = di/a, and F = Xj(T) with probability a^/a. Our second observation is that 
conditional on {T = Tj}, the distribution for T is the same as its unconditional distribution. So, 
F has the same distribution as the random variable obtained by running a Brownian motion X 
with initial distribution P{X(0) = x^} = a,i/a,i = 1, . . . , m, and stopping it independently at 
time T. As a result, F has the desired distribution. 

By conditioning on Z e and letting e — > 0+, the conclusion in the proposition also follows for 
any general initial measure Zq. 

□ 

Remark 4.10. This near extinction behavior is the same as the super Brownian motion (see 
Theorem 1 in [TH]). 

5. Connections with the Arratia flow 

Arratia flow is a stochastic flow which describes the evolution of a continuous family of 
coalescing Brownian motions on R. We refer to for a detailed account and [2] for a survey 
on stochastic flows. By definition, the Arratia flow is a collection {0(s, t,x) : < s < t, x £ M} 
of random variables such that 

• the random map (s,t,x) i— > <f)(s,t,x) is jointly measurable, 

• for each s and x, the map 1 1— > (j)(s,t,x), t > s, is continuous, 

• for each s and f with s < t, the map x i— > <^>(s, t, a;) is non-decreasing and right-continuous, 

• for s < t < u, 4>(t, u, •) o (f)(s, t, •) = (p(s, u, •), 

• for u > 0, (s, t, x) i— > 4>{s + u,t + u,x) has the same distribution as cj), 

• for xi < . . . < x m the process (</>(0, t, xi), . . . , 0(0, t, x m ))t>o has the same distribution 
as a coalescing Brownian motion starting at (xi, . . . , x m ). 

Fix t > 0, it is known that {0(0, i,x) : x G M}, the image of R under map (f)(0,t,.), is a 
discrete set (see pQ). Let ... < x*_ 1 < Xq < x^ < ... be a sequence of random variables such 
that 

(5.1) {0(0, t, x) :xGR} = {..., xl l5 x$, x* l5 . . .}. 

Since Brownian motion has continuous sample paths, the Arratia flow is order-preserving; i.e. 
0(0, t, xi) < . . . < 0(0, t, x m ) whenever xi x m . Set 

IT := sup{x : 0(0, t, x) = x*}. 

Write <f)~ l {^,t,x) for the pre-image of x under 0(0, t, .) Then (Ilj) determines a partition on R 
such that 0" 1 (O,t,x*) = [IT_i,IT[. 
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Not surprisingly, the Arratia flow is closely connected to the process Z studied in the previous 
sections. We first consider its support St- Since St is a discrete set, we can identify it with a 
simple point process by placing a unit mass on each point of St- For any y\ < y 2 < . . . < y 2n , 
by Proposition 14.11 



(5.2) 



{Z t (U r ; =1 ]y 2 ^ 1 ,y 2 ,]) = 0} 



cxp 



jt 



Z (Uj =1 ]Y 2j ^(t),Y 2j (t)]) 



We thus get the following characterization of the avoidance function for St- 

2 



(5.3) P{5 t nU^ 1 j|/ 2i _ 1) |/ 2i ] = 0} 



exp 



it 



Z (U? =1 ]3%_i(t),Yfc-(t)]) 



Consequently the distribution of St is uniquely determined by l|5.3j) : see Theorem 3.3 in |S]. 

(|5.3f) suggests a connection between St and the Arratia flow. Let M t (dy) be a random measure 
on R such that 

2 n 

M t (u" =1 ]y 2 i-i,2/2j]) = — 2 z o(]^(O,t,i/2j-i),0(O,t,y2j)]), yi < 2/2 < • • • < y 2 «- 

Then St can be identified with a Cox process with a finite random intensity measure M t . 
Q5.2JI also leads to a result on the occupation time for Z. For any Borel set 5 in R, 



/ dsP{Z s (£) = 0} = / dsF[exp{-M t (B)}} 
Jo Jo 



A particle representation for Z t is available by using the image of the Arratia flow as a 
skeleton. Given (x*) as in (J5.1|) . let (• • • , £o ; £ii • • •) be independent non-negative random 
variables such that 



'[exp{-A&}| (x*)] = exp 



2AZ o (0- 1 (O,t,x?)) 



2 + A7t 



Then 
(5.4) 

To see this, define 



Zt 



D 



m2" 



where I Q 



m2" 



i=-m2 m 



^ m := ^ m (*)^(0,t,i/2"»), 

i=-m2 m 

is a sequence of independent Feller's branching processes with initial values 



(Z ([(*-l)/2"V/2™[))[' 



m2 m 

i=-m2 m ' 



and in addition, ^"^J is independent of {(p(s,t,x)}. 
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For any aj > 0, j = 1, . . . , n, and y\ < . . . < i/2n, by the same argument as in the proof for 
Theorem 13.21 we have 



lim P 

m— >oo 



lim 

m— >oo 



exp { ~ E a i z t m \]y2j~i,y2-j\) 

3=1 
m2 m 

1 [ exp 



2Z ([(z-l)/2 m ,V2 m [)ELi« J l{]^-i(0,^(0]}(^ 



,!=-m2 m 



lim ] 

m— >oo 



exp 



exp 



exp < 



m2 m 

E 

i=-m2 m 



Z (dx) 



2 + 7 tJ2] =1 a J l{]Y 2j _ 1 (t),Y 2j (t)]}( l /2^) 

2Z ([(i - l)/2 m ,»/2 m [)E; =1 a 3 l{]y 2H (t),r 23 (t)]}(^ 
2 + 7* Ei=i «3 Ul^i-i (t) , ^ (*)]} ( V2 m ) 
2E-=i%l{]^-i(t),^(t)]}(x) 



2 + 7tELi%l{]^-l(t),^2i(t)]}(x) 



^ajZtQyaj-i,^]) 
3=1 



Therefore, 



Further, by the definition of (x*) and the additive property for Feller's branching processes 
we obtain that 



,(m) 



E E trm - E &v 

j=-oo ni_i<j/2 m <rii i=-oo 

Putting these together gives (|5.4j) . 

This interplay is remarkable. On one hand, Z can be constructed using the Arratia flow; on 
the other hand, Z tells us how an initial measure Zq is transported over time under both the 
Arratia flow and the branching. 

The Laplace functional for Zt can also be expressed in terms of (x*) and (H). Given any 
nonnegative bounded continuous function /, for yj = j/2 n , Theorem 13 . 21 yields 

n2 n 



(5.5) 



exp i E f(yj) z t(]yj-i,yj]) 

j=-n2 n 

2 E"=- n2 n /(%)1{]0(O, t, 0(0, t, yj )]}( X ) 



exp 



Z (dx) 



2 + It E"=_„2« /(l/i)W(0, Ifc-i), 0(0, t, yj ))Kx) J _ 



Again, let m — > oo. It follows that 

P exp | -A J f{x)Z t {dx)\ 



exp 



Zo(dx) 



2AE,=-oo/(n,)i{]<,< +1 ]}(x) 
2 + A 7 tE£-oo/(n l )i{]<,< +1 ]}(x) 
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6. A MORE GENERAL MODEL 

Evans observes that what is really at work in the proof for Theorem 13,21 is the additivity 
for the Feller's branching processes. He then suggested that we could use the square of Bessel 
processes (BESQ) to describe the evolution of masses. We are going to carry it out in this 
section. 

For x > and 5 > the square of 6 -dimensional Bessel process starting at x, denoted by 
BESQ 5 (re), is a non-negative valued process £ which solves the following stochastic differential 
equation 

& = x + 2 I ^fisdB s + 6t, 
Jo 

where B is a one-dimensional Brownian motion. The Laplace transform for £ is given by 

(6-1) P [ex P {-A&}] = l — T exp {-y^—} . 

(l + 2At)3 I l + 2Xt) 

Notice that the Feller's branching process is just a BESQ . We refer to Chapter XI in for 
a more detailed introduction on the Bessel processes. 

It is easy to see from (j6.1j) that BESQ <5 (x) is additive in both 5 and x; i.e. if = 1, . . . , m} 
is a sequence of independent processes such that each £j is a BESQ*(xi). Then YZi& 

is a 

BESQ&^QXiZi)- 

Now we are going to modify the process Z defined in Section 3 by letting the masses of the 
particles be governed by the BESQ processes. Since the dimension is an additional parameter 
for BESQ, we need to introduce another measure- valued process A to describe the evolution of 
the dimension. 

As in Section 3, we first consider two systems of interacting particles. Given a finite measure 
Zq on R, for any m, choose x\, . . . , x m to be i.i.d. random variables with a common distribution 
Zq := Zq/Zq(M). Let (X%, . . . ,X m ) be an m-dimensional coalescing Brownian motion starting 
at (xi , . . . , x m ) . 

Given another finite measure Ao on M + with a finite "mean" [i := <5Ao(ci<5), let 6%, . . . , 5 m 
be i.i.d. random variables with a common distribution A := Ao/Ao(M). We further suppose 
that and (Aj) and (5i) are independent. 

Put 5 := A (R) and z := Z (R). Let (&,..., £ m ) be a collection of m independent BESQ^/ m (z/ 
processes. Then 

m 
i=l 

and 

i=l 

define two M^(R)-valued processes. 
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Similar to Lemma 13, II we can show that both {Z( m ^} and {A( m )} are C-relatively compact in 
D(Mp(M)). They have unique weak limits by Theorem 16.11 which we will prove shortly. 

Let Z and A be the weak limits for {Z^} and {A( m )}. Intuitively, {(Z (B), A (B)) : 
B G SQR)} describes the initial mass-dimension distribution on R, and {(Zt(B), At(B)) : B £ 
i3(]R),0 <t< 00} describes the simultaneous mass-dimension evolution for such a model, which 
we call a super square of Bessel process with spatial coalescing Brownian motion. 

For any nonnegative constants ctj,(3j, j = 1, . . . , n and t > 0, put 



I t (x) :=Y,^{\Y2^i(t),Y 2] (t)]}^ 



and 



J t {x) :=Y^m]Y2 3 -i{t),Y 2j (t)}}{x), 

where, as usual, (Yi, . . . , Y-2 n ) is an 2n-dimensional coalescing Brownian motion starting at 
{yxi ■ ■ ■ 1 U2n)- The next result determines the joint distribution for (Zt(B),At(B)), B G 

Theorem 6.1. For any ctj > 0, > 0, j = 1, . . . , n and any yi < . . . < y 2n , we have 



(6.2) 



exp ^ -^2a j Z t (]y 2 j-i,y2j}) - ^ (3jA t (]y 2j -i, V2j]) 
3=1 3=1 



exp 



Zo(dx) 



£ In (1 + 2tl t (x)) + - /ff + ±J t (z) 
2z l + 2tl t (x) z 



Proof. To prove (|6.2|) . again, we first fix ... , 5 m ) and (£1, . . . , Cm)- It follows that 
(6.3) 

n n 



3=1 



3=1 



n 1 m 



1=1 



3=1 \i=l '" 



n 1 m 



ex P {-Y,[zZ a Mt)H]Y2 j -i(t),Y 2j (t)]}(x l ) 

3=1 \i=l 

• £-^l02 ; .!/). V 2/ !7);} ( , ! 



1=1 



m 



6i5 



i=l 



i=l 
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We then fix (x{) and (Si), and take expectations with respect to (£1, . . . , £ m ). By (|6.1|) the right 
hand side of ()6.3j) is equal to 

(6.4) 



Y[(l + 2tl t (xi)) 2™ exp 



i=l 



zl t (xi) 



SiSJ t (xi) 



m(l + 2tI t ( Xi )) 



m 



Z (dx) / A (dS)(l + 2tI t (x))~^ exp 
Jo 

Let m — > oo in (|6.4j) . We finally have 

f n n 

exp < -^a j Z t (]y 2 j-i,y2j]) - ^ /3j A t (]y 2 j-i, V2j\) 



zL(x) 



SSJ t (x) 



m(l + 2tl t (x)) 



m 



3=1 



3=1 



lim P 

m— >oo 



lim P 

m— >oo 



lim P 

m— >oo 



j=l 3=1 

ss 



exp 



exp 



Ixl-l 



— oo 
oo 



Z (dx)Ao(cW) ( 1 - — In (1 + 2tl t (x)) ) [ 1 



SS 



zl t (x) 



SSJ t (x) 

m(l + 2tlf(x)) m 
zl t (x) SSJ t (x) 



Z ( dx) I M d0 I 1 - -ln(l + »/,(,)) - + ^ 
Z„(dx) / A„(dJ) I y In (1 + 2(/,(:c)) + ; + fflV,^ 



y°° Z (dx) ^ In (1 + 2tJ t (x)) + - 



h(x) 



+ 2tl t (x) z 



□ 



Remark 6.2. Notice that A is just the process Z in Theorem 13 . 2 1 wit h 7 = 0. 



The generalized model considered in this section will not die out for /i > 0. Many of the 
properties in Section 3 and 4 can be discussed in a similar fashion. But we leave the details to 
the interested readers. 

Acknowledgement: The author is grateful to Steven Evans for a suggestion that results in 
Section 5 of this paper. The author also thanks Carl Miiller for a helpful comment. 
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